Noise-driven transition to quasi-deterministic limit cycle dynamics in excitable 

systems 
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The effect of small-amplitude noise on excitable systems with large time-scale separation is ana- 
lyzed. It is found that small random perturbations of the fast excitatory variable result in the onset 
of a quasi-deterministic limit cycle behavior, absent without noise. The limit cycle is established 
at a critical value of the amplitude of the noise, and its period is nontrivially determined by the 
relationship between the noise amplitude and the time scale ratio. It is argued that this effect might 
provide a mechanism by which the function of biological systems operating in noisy environments 
can be robustly controlled by the level of the noise. 
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Understanding the effect of random perturbations in 
dynamical systems is a fundamental problem with a large 
number of applications in physical, chemical, and biolog- 
ical sciences. When these perturbations are small, noise- 
driven systems may often exhibit rare events which lead 
to switching- like dynamics on long time-scales This 
behavior is well understood in systems obeying detailed 
balance, which is typical of systems close to thermal equi- 
librium (see, for example, [2j). 

Systems far from equilibrium often fail to obey de- 
tailed balance. Even in the absence of noise, these sys- 
tems are capable of having complex dynamical behav- 
iors. One example of such a behavior is limit cycle os- 
cillations. Another example is excitability: in excitable 
systems, small perturbations of the unique stable steady 
state decay, while sufficiently large perturbations lead to 
large-amplitude dynamical excursions before the system 
returns to the steady state (see, for example, 0). Ex- 
citable systems are very common in biology, with nerve 
cells being just one example (see, for example, Q). 

The purpose of this Letter is to show that under cer- 
tain conditions small noise may transform the dynamics 
of an excitable system into a quasi- deterministic limit 
cycle, signifying a transition from excitable to oscillatory 
behavior under the action of the noise. This transition 
occurs at a critical value of the noise amplitude. Further- 
more, the frequency and other parameters of the limit 
cycle are controlled by the amplitude of the noise, which 
distinguishes this mechanism from other mechanisms of 
noise-induced coherence (see, e.g., [HEIB])- 

The key to our argument is the existence of a strong 
time scale separation in the system's dynamics. In biol- 
ogy, this is typically the case: excitability often arises as 
a result of the competition of positive and negative feed- 
backs operating on different time-scales, with fast exci- 
tatory and slow recovery variables 0]. In this situation 
it is possible to have an interplay between the slow de- 
terministic time-scale of the recovery variable and the 



exponentially long time-scale associated with rare events 
induced by the noise. It is precisely this interplay that 
makes the considered effect possible. 
Consider a general dynamical system 



v = ag(u,v). 



(1) 
(2) 



In the context of excitable systems, u and v (generally, 
vectors) will be sets of excitatory and recovery variables, 
respectively; / and g are the nonlinearities, a is the ratio 
of the time-scales, rj is some external noise perturbing the 
excitatory variables, and e <C 1 measures its amplitude. 
If a <C 1, there is a large time-scale separation in the 
deterministic part of the dynamics governed by / and g, 
with u and v being the "fast" and the "slow" variables, 
respectively ^jj- That is, on the time-scale of order I the 
dynamics of u is governed by Q with v frozen. There- 
fore, in the absence of the noise, the trajectory quickly 
approaches the slow manifold defined by the solutions of 
f(u, v) = for fixed v, and then proceeds on the time 
scale a -1 along this slow manifold and into the globally 
attracting equilibrium point. 

With the introduction of the noise, the situation 
changes. Indeed, the trajectory may escape the slow 
manifold in the u direction via a noise-activated process 
on the Arrhenius time-scale t = v^ 1 exp(j3(v)/e) ^S> 1, 
where /3(v) is some energy barrier to be crossed to initi- 
ate the escape and v is some characteristic frequency Q . 
As a result, the system can perform a large excursion 
driven by the deterministic part of the dynamics, after 
which the trajectory can land again somewhere on the 
slow manifold. The dynamics can then proceed along 
the slow manifold until another escape event happens, 
and so on. This process may lead to a bona-fide limit 
cycle because of the following two ingredients. First, the 
interplay between the escape events and the motion along 
the slow manifold requires that their time-scales be com- 
parable. But since the escape rate is a rapidly varying 
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FIG. 1: Deterministic flow generated by the nonlinearities 
from Eq. (QJ. Results of the numerical solution of Eqs. Q 
and in the u — v plane for A = 0.7 and a — 0.01. Thick 
red and blue curves are the u- and ii-nullclines, respectively. 

function of v, escape from the slow manifold will occur 
with overwhelming probability in the small vicinity of 
point v — v* on this manifold, where satisfies 

/3(0=e log a~\ (3) 

Here we assumed that a,e — > in a way that elogcv -1 
remains finite and neglected higher order terms. The 
next ingredient to obtain a limit cycle behavior is then a 
mechanism of reset. It requires that after a large excur- 
sion initiated from v+, the trajectory returns to the slow 
manifold to a point which leads again to by the slow 
motion. Then the process will repeat itself indefinitely 
and the dynamics of the system will indeed be a limit 
cycle. Let us point out that the mechanism of escape 
at a fixed v+ is related to the phenomenon of stochastic 
resonance (for a review, see Q). Also note that the pos- 
sibility of limit cycle behavior under the action of small 
noise in systems lacking such behavior in the absence of 
the noise was recently pointed out by Freidlin ||. 

To demonstrate the feasibility of the scenario above, 
we will consider the Brussclator This is a prototyp- 
ical model of a system far from equilibrium in which the 
nonlinearities take on the following form 

/ = 1 + u 2 v - (1 + A)u, g = Au- u 2 v. (4) 

Here u and v are scalars and A is a control parameter. 
For simplicity we will assume that r\ is external white 
Gaussian noise, with (rj(t)r)(t')) = S(t — t'). 

The Brusselator is an excitable system when a -C 1 
and A < 1. This can be seen from its phase portrait 
shown in Fig. ^ When A < 1, the nullclines of Eqs. 
and (0) intersect on the stable branch of the it-nullclinc, 
so the flow is always into the unique equilibrium point 

uq = 1, vq = A. (5) 

Note that the slow manifold here is the rising part of the 
M-nullcline (see Fig. . It is also clear from the figure 



that sufficiently large increases in the u variable away 
from equilibrium will result in large excursions. 

We performed numerical simulations of the stochas- 
tic differential equations given by Eqs. QJ, J3J), and 
Fig. |2Ja) shows the time series for one realization of the 
noise. This figure shows a train of large amplitude spikes 
in the excitatory variable. What is striking is that the 
spikes are occurring in an almost periodic fashion, with 
their amplitude and other characteristics being approxi- 
mately the same. This can be better seen from the phase 
plane, Fig. EJb). In fact, Fig. 0b) strongly suggests a 
nearly limit cycle behavior in the u and v variables. To 
verify this, we computed the stationary probability dis- 
tribution function (PDF) for this process (not shown). 
The PDF is found to be concentrated around a part of 
the slow manifold. One can also see a pronounced loop in 
the PDF for large values of u. Importantly, the PDF has 
a maximum at some value of v = v* < vq- Let us point 
out that this result shows that our system lies beyond 
the region of applicability of the large deviations theory 
(if the latter were true, the PDF would be monotonically 
decreasing away from the equilibrium point 0). 

To further investigate the effect of the noise, we have 
analyzed the statistics of the interspike time intervals for 
different values of the noise amplitude e. For the purpose 
of this analysis, we defined as a spike any excursion with 
an amplitude u max > 10. Fig. |21 shows the mean inter- 
spike distance T and its standard deviation ctt obtained 
from the numerical solution of Eqs. JIJ, 0) and 0) for 
different value of the amplitude of the noise. Also, in the 
inset we show the ratio <jt/T as a function of e, which 
characterizes the "signal-to-noise ratio" for the interspike 
distance. 

First, observe that for very small e the spikes are ex- 
ceedingly rare and have the character of a Poisson pro- 
cess (see the inset). They represent rare large-amplitude 
fluctuations away from the equilibrium point. Then, at 
some critical value of e ~ 0.02 the increase in the noise 
amplitude results in a rapid decrease in the interspike 
distance. Furthermore, the ratio <jt/T rapidly decreases 
and stays low in a broad range 0.05 < e < 0.5, signifying 
high degree of signal coherence (see Fig. EJ) . Let us em- 
phasize that in this range T shows significant dependence 
on e, while ctt does not (as seen from the errorbars). In 
other words, the noise amplitude really acts as a control 
parameter for the deterministic oscillatory behavior of 
the system. Finally, for larger noise amplitude, the spike 
train looses coherence again, since then the noise is no 
longer weak. 

To explain these observations and corroborate the gen- 
eral scenario given above, let us consider the limit of 
strong time-scale separation. On the fast time-scale the 
value of v is fixed, hence Eq. (JJJ describes the motion 
of a particle in the potential well of the form V(u) = 
— ^vu 3 + |(1 + A)u 2 — u. This is the classical escape 
problem of Kramers, for which the average escape time 
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FIG. 2: Numerical solution of the stochastic differential equation. The parameters are: A = 0.7, a = 0.01, e = O.f. (a) The 
time series with the quasi-periodic spike train, (b) The phase plane plot showing the quasi-deterministic limit cycle. Red and 
blue lines are the it- and u-nullclines, respectively. Solid green line is the predicted escape path and dashed green line is the 
excursion following the escape. 



is given by the following formula 0] 

2» na + Af 
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(6) 



From this equation, we obtain explicitly (3(v) — 
[(1 + A) 2 - Av] 3 / 2 /3v 2 . This f3{v) is inserted into Eq. ©, 
which fixes function of a and e. 

After the trajectory escapes the neighborhood of the 
stable nullcline at v = u*, it continues moving toward 
increasing values of u. At this point the effect of the 
noise becomes negligible. With the increase of u, the 
effective time-scale of v decreases [see Eq. (j2J], so when 
the system undergoes a large excursion, the time-scales 
of these two variables can no longer be separated. On 
the other hand, when u ^> 1 and v ~ 1, Eqs. Q and (J2J 
can be simplified by neglecting all the terms except u 2 v. 
The resulting system of equations with the asymptotic 
boundary conditions u{— oo) = 0, v{— oo) = can be 
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FIG. 3: Mean interspike distance T and its standard deviation 
<jt (shown as errorbars) as a function of the noise amplitude. 
Inset: the plot of the ratio of the standard deviation to the 
mean. Except for the noise, the parameters are as in Fig. 



integrated to give v = v+ — au and it — it 2 (v* — cm). 
This equation can be solved exactly to give the transition 
layer for the rising part of the spike. It shows that in the 
spike u rises to u max = a w* > 1 on the time-scale of 
a <C 1, while v approaches zero asymptotically. This is 
then followed by the fall of the trajectory onto the u- 
nullcline with fixed (small) v on the slower time scale of 
order 1, see Eq. with v — (a more precise layer 
analysis of the solution can be performed [111]). 

Following the excursion, on the slow time-scale the tra- 
jectory moves along the u-nullcline, starting at v = 
(asymptotically) . After a little algebra we obtain for the 
slow motion: 



2a(A 



1 + A - 2v + + A) 2 - Av 



(7) 



From Eq. (JJJ, v > for < v < vq, so the trajectory will 
reach the point v — on the u-nullcline in time T = 
T(v*), which can be explicitly calculated by integrating 
Eq. (JJJ. Upon reaching v*, the trajectory jumps off from 
the nullcline and returns onto the u-nullcline at the point 
where v = 0, thus completing the loop [cf. Fig. Hfb)]. In 
summary, we obtain a limit cycle whose amplitude and 
period are asymptotically determined by the value of v+, 
which is in turn a function of a and e. 

Since the system spends most of the time on the slow 
manifold, asymptotically the period of the limit cycle will 
be equal to T(u*). This prediction can be compared with 
the mean interspike distance obtained from the simula- 
tions of Fig. We found that the two are in qualitative 
agreement for 0.05 < e < 0.5, but T(v*) consistently un- 
derestimates the observed value of T by a factor of about 
1.5. This is not unexpected, since the small parameter 
of the asymptotics is l/lna^ 1 , which in practice is not 
very small. We verified that the accuracy increases with 
decrease of a. 

The above analysis also predicts the existence of a crit- 
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ical amplitude of the noise for the establishment of the 
limit cycle behavior. Indeed, the attainable values of v 
on the u-nullcline lie in the interval < v < v . There- 
fore, the solution of Eq. J3J may lie in this interval only 
if e > e r , where 



Er. = 



3A 2 log cr 



(8) 



No limit cycle behavior is possible when e < e c . As 
e approaches e c from above, we have u* — > vo, and 
T(v*) — > oo. Therefore, for fixed deterministic part of 
the dynamics there is a transition to a limit cycle be- 
havior at a critical value of the amplitude of the noise. 
We verified that the predicted value of e c gives the cor- 
rect order of magnitude for the onset of the oscillatory 
behavior. 

Let us emphasize that the mechanism described above 
is robust: it does not require fine-tuning of the system's 
parameters. This is in contrast with other mechanisms, 
by which noise can induce coherence (see, for example, 
Those mechanisms require that in the absence 
of the noise the system is near the threshold between 
excitable and oscillatory behavior. 

Now we discuss the potential implications of the ob- 
served phenomena. First of all, our results suggest that 
in systems with strong time-scale separation one should 
be careful in identifying the class of dynamical models 
for explaining the observations of oscillatory behaviors. 
Real systems are always noisy, and our analysis indi- 
cates that oscillatory behavior can be obtained in systems 
with intrinsically non-oscillatory dynamics. In particu- 
lar, this may be relevant to the identification of cell types 
in model neural networks (see for example the discussions 
in [l2l Il3jp. We note that the noise-induced transition 
similar to the one discussed in this Letter was recently 
observed in the numerical simulations of randomly forced 
Hodgkin-Huxley neuron [T^j . 

Another possible implication has to do with coupled 
excitable systems. Our analysis indicates that in such 
systems the level of the noise, both extrinsic and intrin- 
sic, may be used as an information carrier and trans- 
formed into (quasi-)deterministic signal. As a prototype, 
consider a system of all-to-all positively coupled excitable 
cells. Under the action of the noise of sufficiently small 
amplitude each cell will occasionally generate a spike. 
These spikes will have random phases, so their total in- 
put on each individual cell may average to a stationary 
random signal of low intensity. Now, if the noise level 
suddenly increases due to an external disturbance, the 
cells may switch to the noise-assisted oscillatory mode. 
This will further increase the effective noise amplitude, 
so that the oscillatory mode may persist even after the 
disturbance is removed. In colloquial terms, the system 
in a dormant state may wake up from the outside rattle. 

In a similar way, our results may be applied to spatially 



distributed excitable media [3j . In these systems the ana- 
logue of the noise-activated event will be the formation 
of the radially-symmetric nucleus, leading to the conse- 
quent initiation of the radially-divergent wave 

[303. in 

the wake of such a wave the system will undergo recovery. 
It is clear, then, that the system will be most recovered 
at the position where the wave was initiated. Hence, the 
new wave will be initiated again at the same spot, with 
the dynamics repeating periodically. This suggests that 
the well-known phenomenon of target pattern formation 
in two-dimensional excitable media Q might have an al- 
ternative interpretation in terms of noise-driven quasi- 
periodic wave generation. 
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